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I.  BACKGROUND 


S 


This  work  was  supported  by  the  Office  of  Naval  Research  under 
Contract  N00014-73-C-0023,  ONR  Contract  Authority  Identification 
Number  NR  384-320/04-27-72/468. 

II.  DISCUSSION 

The  work  reported  here  is  a  continuation  of  the  research  con¬ 
ducted  under  our  former  Contract  N00014-71  -C-0050.  Both  theoretical 
and  experimental  work  has  been  undertaken.  Progress  in  ceveral  areas 
is  summarized  below. 

A.  Attenuation  Compensation  for  Measuring  Ultrasonic  Velocity 
at  High  Temperatures  in  One  or  More  Segments  of  Wire 

Since  attenuation  increases  with  temperature,  particularly  above 
half  the  melting  point  (absolute  scale),  it  is  necessary  to  compensate  for 
attenuation  to  assure  echoes  of  large  enough  amplitude  to  measure  when 
working  with  wires  immersed  in  a  hot  environment.  A  method  iB  out¬ 
lined  to  calculate  the  echo  amplitudes  desirable  at  room  temperature  in 
order  to  fabricate  a  wire-sonic  transmission  line  for  operation  at  high 
temperatures  for  velocity  measurements.  Applications  to  ultrasonic 
thermometry  are  mentioned. 

B.  Versatile  Ultrasonic  Facility  for  Operation  to  1000°C 

To  facilitate  measurements  on  wires  and  thin  rode  from  room 
temperature  to  1000  C,  a  simple  and  versatile  vacuum  furnace  and  heat 
controller  have  been  built.  An  ultrasonic  port  allows  the  entrance  of 
several  wire  transmission  lines  and  a  thermocouple.  Ultrasonic  mea¬ 
surements  of  modulus  and  attenuation  can  be  made. 

C.  Measurements  on  Remendur  to  1000°C 


In  the  vacuum  furnace  mentioned  above,  measurements  of  Young’s 
modulus  and  the  attenuation  of  extensional  ultrasonic  waves  have  been 
made  to  1000  C  in  Remendur,  a  magnetostrictive  iron/cobalt  alloy.  It 
was  desired  to  study  magnetostrictive  metals  up  to  and  beyond  the  Curie 
temperature.  Further  measurements  are  planned. 
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D.  Statistical  Errors  in  Attenuation  and  Reflection  Coefficient 

Statistical  methods  have  been  used  to  determine  the  errors  to  be 
expected  in  measuring  the  attenuation  a  in  a  specimen  of  length  L  and 
the  reflection  coefficient  R  at  a  buffer/ specimen  interface  (or  lead-in¬ 
wire/wire -specimen  joint)  with  traveling  wave  methods.  Formulas  are 
derived  for  the  maximum  (3a)  errors  in  R  and  2orL. 

E.  Momentary  Contact  Experiments  for  Moduli  and 

Attenuation  at  Elevated  Temperatures _ 

We  are  preparing  to  experiment  further  with  superalloys  at  elevated 
temperatures.  The  next  alloy  to  be  tested  will  be  Waspalloy,  a  structural 
metal  used  at  high  temperatures. 

F.  Polymerization  of  Plastics 

) 

Our  work  on  testing  the  curing  of  epoxy  resins  has  been  written 
up  for  technical  publication;  the  text  appears  as  Appendix  A.  The  paper 
has  application  to  various  plastics  and  coatings,  and  possibly  to  solids 
during  melting  and  freezing.' 

G.  Nomogram  for  Attenuation  and  Reflection  Coefficient 

To  familiarize  the  scientific  and  engineering  community  with  the 
utility  of  making  ultrasonic  attenuation  measurements  by  ultrasonic 
traveling  wave  methods  using  buffer  rods  or  lead-in  wires,  a  letter  has 
been  prepared  for  publication  showing  a  nomogram  for  obtaining  the  at¬ 
tenuation  and  the  reflection  coefficient  from  the  amplitudes  of  three  echoes. 
This  letter  is  given  in  Appendix  B. 

H.  Ultrasonic  Velocity  in  Sheet  Material  ' 

Ultrasonic  velocity  in  sheet  material  can  be  used  to  measure  or 
monitor  the  moduli  of  the  material.  Various  production  variables  might 
be  brought  under  closer  control  if  the  moduli  of  sheets  could  be  monitored. 
A  set  of  experiments  was  performed  on  various  types  of  paper  to  investi¬ 
gate  wave  propagation  in  a  material  both  high  in  attenuation  and  highly 
anisotropic.  A  paper  was  written  summarizing  these  measurements  and 
reviewing  some  foreign  and  domestic  work.  This  paper  appears  as 
Appendix  C.  It  is  f&it  that  other  types  of  sheet  material  could  be  monitored 
at  high  temperatures  by  ultrasonics. 
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HI.  DETAILS  OF  WORK 

A.  Attenuation  Compensation  for  Measuring  Ultrasonic  Velocity 
at  High  Temperatures  in  One  or  Mora  Segments  of  Wire 

In  certain  experiments,  one  may  be  interested  in  measuring  the 
ultrasonic  propagation  constants  along  one  or  more  segments  of  a  wire 
waveguide.  For  example,  one  may  wish  to  compare  the  response  of  two 
or  more  different  materials,  e.  g.  ,  pure  Re.  pure  Mo,  and  Re/Mo  alloys, 
in  a  given  thermal  and/or  nuclear  environment.  Alternatively,  one  may 
wish  to  use  a  given  material,  e.  g.  ,  Re,  to  sense  the  variation  or  profile 
of  an  environmental  parameter  such  as  temperature  or  neutron  fluence. 

At  temperature  below  about  half  the  absolute  melting  temperature, 

T  <  0.  5  T^,  the  attenuation  coefficient  a  at  frequencies  well  below  1  MHz 
is  usually  small  enough  so  that  detectable  echoes  can  be  obtained*  even 
when  the  sound  pressure  reflection  coefficient  R  is  much  less  than  unity; 
e.  g.  ,  R  ®  0. 1.  In  this  case  at  least  10  and  possibly  20  segments  can  bv 
monitored  by  a  single  transducer. 

On  the  other  hand,  at  temperature  above  ~0.  5  Tj^,  a  may  increase 
by  more  than  1  order  of  magnitude.  To  provide  sufficiently  strong  echoes, 
one  may  have  to  divide  the  available  ultrasonic  energy  among  a  smaller 
number  of  segments,  e.  g.  ,  5  or  less.  Correspondingly,  the  R's  at  seg¬ 
ment  interfaces  need  to  be  increased. 

The  problem,  then,  is  to  optimize  the  choices  of  Rj,  for  a  given 
set  of  a  L-*s,  where  the  Lj's  are  the  lengths  of  the  segments. 

Let  us  define  "optimum"  to  be  that  set  of  Rj  and  L^  which  yields 
essentially  equal  magnitude  echoes  at  the  expected  conditions  of  maximum 
attenuation,  and  minimum  or  tolerable  interference  echoes  at  temperatures 
well  below  0.  5  Tj^j.  To  illustrate  one  approach  to  this  optimization  pro¬ 
cedure,  let  us  analyze  the  case  where  temperature  is  uniform  over  the 
segments,  and  high  enough  above  0.  5  T^  so  that  very  small,  identical 
R's  such  as  R  =  0.  1  will  not  suffice.  Further,  since  moderately  large 
R’s  will  be  large  enough  to  yield  interfering  reverberations,  the  lengths 
of  adjacent  segments  need  to  be  unequal.^'  * 

The  temperature  dependence  of  or  in  Re  is  shown  in  Fig.  1.  It  is 
seen  that  a  increases  sharply  above  about  1600°C  from  0.  005  Np/cm 
(4.  34  dB/m)  to  0.  03  Np/cm  (26  dB/m)  at  2000°C.  An  a  of  0.  015  Np/cm 
(10  dB/m)  occurs  near  l650°C.  On  a  normalized  (homologous)  tempera¬ 
ture  scale,  l650°C  translates  to  *  0.  56  Tj^j. 
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In  Fig.  2  we  plot  lines  of  positive  slope,  for  several  a* s.  Let 
us  use  this  graph  to  determine  R^.  For  simplicity,  suppose  we  are  to 
operate  at  a  maximum  temperature  of  l650°C,  i.  e.  ,  a  =10  dB/cm. 

Further,  suppose  we  are  to  utilize  a  waveguide  containing  four  segments 
of  alternating  lengths  =  30,  20,  30,  20  cm,  as  shown  in  Fig.  3.  If 
the  echoes  A,  B,  C,  D,  E  are  to  be  equal  at  l650°C,  their  relative 
amplitudes  at  room  temperature  may  be  obtained^frgm  Fig.  2,  remember¬ 
ing  that  each  segment  attenuates  the  signal  by  e  1  :  0. 199,  0.  316,  0.  631, 
and  1.  000  V. 


Starting  with  echo  E  at  the  free  end,  we  may  proceed  to  compute 
Rj  and  then  compute  corresponding  diameter  ratios. 

Results  are  tabulated  in  Table  I.  A  comparison  with  cases  where 
a  is  negligibly  small,  as  well  as  where  a  is  substantially  larger  than 
10  dB/m,  is  given  in  Table  II. 

In  principle,  a  rod  containing  diameter  steps  as  required  by  Table  I 
could  be  fabricated  by  centerless  "plunge"  grinding;  by  welding  together 
segments  of  the  required  diameters;  by  electric  discharge  machining;  by 
drawing;  or  by  other  means.  Figure  4  is  included  to  illustrate  that 
stepped- diameter  wires  can  in  fact  be  built,  and  when  heated,  provide 
a  detectable  echo  pattern.  However,  once  we  have  calculated  the  Rj’s, 
we  may  produce  the  reflections  by  means  other  than  stepping  the  diameters 
of  each  complete  segment.  For  example,  circumferential  grooves  or 
added  masses  such  as  weld-upset  bumps  may  be  used  (Fig.  5).*  In  some 
materials,  another  alternative  is  drilling  axial  holes  of  appropriate 
diameters  and  depths. 


Details  on  the  calculation  of  R-  and  diameter  ratios  are  given  next, 
using  the  notation  of  Fig.  3. 

t 

For  the  end  segment, 

TJR* 

D  _ D  1  D  0.631  ,  . 

E  '  RE  T34T43  '  T34T43  =  10 


since  Rg  = 
But 


-  1  at  the  free  end. 


r34-  1 
r34+1 


(2) 
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Parameters  calculated  for  segmented  line  of  one  material 


Segment 

1  2  3 

Diameter 

Notation  D 

D1  °2  D3 

Relative  (mm)  2.  36 
Diameter  (in.  )  0.  0943 

2.30  2.05  1.69 

0.0918  0.0818  0.0675 

Length  L.  (cm)  - 

30  20  30 

Round-Trip 
Attenuation  (dB)  - 
if  a  =  10  dB/m 

6  4  6 

Required  Relative 
Echo  Amplitude 
at  Room  Temp,  V 

0.  100 

0.199  0.316  o.i 

Sound  Pressure 

Refl.  Coef.  R 

0*  0262 

-0.1175  -0.19  -0.- 

Transmission 

Coef.  T* 

0.973 

0.  883  0.  81  0. 

Transmission 

5 


Attenuation 


Low 


Medium 


High 


Typical  |  R| 

<0.  1 

~0.  1  to  0.  5 

>  0.  5 

rsj 

No.  of  segments 

>  id 

<  5 

2 

Relative  lengths 

Equal 

Adjacent 
segments 
unequal,  but 
alternate 
segments 
may  be  equal 

Unequal 

Mode 

Extensional  or 
torsional 

Extensional  or 
torsional 

Extensional 
and/ or 
torsional 

For  the  next  segment  we  similarly  write 


B 

C 


T  T 
12  21 


0. 199 
0.  316 


0.  63 


(14) 


which  leads  to  D./D  =  0.  891,  R  =  -  0.  1175,  T,  _  =  0.  883,  T._  =  1.  115. 

1  13  l  c  d»  1 


For  the  fourth  segment, 


\ 


A  fA _ 1 _ 0.  100 

B  “  rb  toitio  '  °‘199 


(15) 


from  which  there  results 


D  /D  =  0.975 
1  o 


(16) 


The  Remendur  transducer  diameter  may  be  made  large  enough  to 
impedance  match  to  the  Re  diameter  Dc,  or,  alternatively,  the  Re  may 
be  tapered  down  to  match  a  given  Remendur  diameter. 

5  6 

Details  on  transducer  construction,  self-heating  of  Re  to  T^, 
and  scanning  segmented  lines  to  measure  selected  echo  pairs7  have  been 
described  elsewhere,  and  so  are  not  repeated  here. 

From  Fig.  3  and  Table  I  it  is  seen  that  the  magnitude  of  the  R^'s 
increases  monotonically  from  0  at  the  transducer/lead-in  joint  to  unity 
at  the  free  end.  From  the  tabulated  transmission  and  reflection  coef¬ 
ficients,  and  assuming  a  =  0  and  then  a  =  10  dB/m,  we  can  calculate  the 
amplitudes  of  the  interfering  reverberations  such  as  B’,  C',  5",  C".  The 

-2aL  -Zab^ 

formulas  are  of  the  form  B’  =  -  BR^R^  e  1  and  C'  =  CR^R^.  e 

B*  and  C*  are  inverted,  and  are  straddled  by  C  and  D.  B"  and  C”  are  not 
inverted,  and  under  isothermal  conditions,  coincide  and  therefore  add, 
midway  between  D  and  E.  Their  values  relative  to  E  =  + 1  V  are  given 
in  Table  III.  It  is  seen  that  these  reverberations  are  smaller  than  the 
next  desired  echoes  by  about  2  or  more  orders  of  magnitude.  Therefore 
the  diameter  ratios  calculated  above  satisfy  the  given  optimization  criterion. 
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Table  III 


Regarding  Fig.  5,  we  note  that  when  the  coil  (transducer)  is  near 
the  center  of  the  line  containing  circumferential  "bumps'1  and  "grooves,  " 
(see  bottom  oscillogram)  the  echoes  show  corresponding  polarities  from 
both  sets  of  discontinuities.  However,  when  the  coil  is  moved  to  either 
end,5  the  echoes  from  the  closer  set  may  be  lost  in  the  recovery  transient, 
and  so,  after  recovery,  one  observes  a  set  of  like  echoes  from  the  remote 
segments.  Again,  these  echoes  are  due  to  |  Rj's  small  enough  so  re¬ 
verberations  in  the  segments  are  small  enough  not  to  interfere.  While 
the  oscillogram  traces  in  Fig.  5  were  generated  using  a  Remendur  dem¬ 
onstration  line  of  length  less  than  1  m,  one  could  construct  a  much  longer 
line,  ~10  to  100  m,  if  clearer  separation  of  bipolar  echoes,  detected  with 
a  near-center  transducer,  were  desired. 

We  may  further  note  in  Fig.  5  an  "apparent"  polarity  reversal 
for  bumpy  or  groovy  echoes,  when  the  coil  is  moved  from  either  end  to 
the  center.  This  illusion  is  due  to  what  appears  as  a  phase-reversing 
effect  of  the  reinforcing  end  echo,  when  the  coil  is  at  the  end.^  Careful 
examination  of  the  slope  of  the  leading  edge  of  the  echoes  will  show, 
however,  that  the  leading  edge  retains  its  polarity,  when  the  coil  is 
moved  from  the  center  to  the  end  of  the  line. 


B.  Versatile  Ultrasonic  Facility  for  Operation  to  1000°C 

The  system  for  operating  to  1000°C  consists  of  an  evacuated 
quartz  glass  tube  heated  at  its  central  portion  by  a  clamshell  furnace. 

The  furnace  (Hevi  Duty  Electric  Company  Type  M-2012)  is  controlled 
manually  with  20  Amp  Variac  connected  as  shown  in  Fig.  6.  The  power 
source  was  mounted  in  a  portable  console  for  safety  and  convenience. 
Temperature  is  monitored  with  a  chromel  alumel  thermocouple  con¬ 
nected  to  a  Leeds  &  Northrup  potentiometer.  Inside  the  furnace,  a 
steel  cell  is  built  to  keep  even  temperature  over  the  entire  volume. 

The  quartz  glass  tube  shown  in  Fig.  7  is  evacuated  by  a  roughing 
pump.  Baffles  are  put  on  both  ends  of  the  tube  to  keep  the  radiant  energy 
from  damaging  the  vacuum  seals  at  the  ends  of  the  tube. 

Tungsten  v/ires  are  used  as  lead-in  which  go  through  a  Conax 
vacuum  seed  into  the  cell  as  described  in  Fig.  7.  Inside  the  cell  a  steel 
rack  is  used  to  support  the  samples.  Very  thin  steel  wire  is  used  to 
hold  each  tungsten  lead-in  wire  in  position  with  respect  to  the  rack.  We 
try  to  have  as  small  a  contact  of  the  tungsten  lead-in  with  the  rack  as 
possible  so  that  we  won't  have  any  noise  (spurious  reflections)  caused 
by  the  rack.  The  rack  with  wires  is  described  in  Fig.  8. 

Each  specimen  is  welded  to  a  lead-in  line  of  tungsten  which  is 
terminated  outside  the  glass  tube  by  a  magnetostrictive  Remendur  stub. 
The  lines  are  driven  by  Panametrics'  magnetostrictive  transducer 
Model  KT-55  and  Panametrics'  Panatherm(^y  Model  5010,  and  monitered 
with  a  Tektronix  535  oscilloscope. 

In  the  whole  experiment,  we  keep  the  vacuum  below  90  micron 
in  order  to  minimize  the  oxidation  of  the  samples  and  lead-in. 

C.  Measurements  on  Remendur  to  1 000  C 


1.  Samples 

Several  samples  of  Remendur  wire  have  been  measured.  We  are 
reporting  on  two  of  them. 

The  dimensions  of  Sample  3  and  Sample  4  are  given  in  Table  IV. 
The  lead-in  wires  are  tungsten  of  0.  1  cm  in  diameter. 
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Sample 


Table  IV 
Diameter 


Length 


3 

4 


0.  075  cm 
0.  075  cm 


7.  04  cm 
9.  32  cm 


a5i 


Samples  3  and  4  are  simultaneously  in  the  furnace  through  the 
measurements,  Run  3  up,  Run  4  up,  and  Run  4  down,  for  longitudinal 
modulus  and  attenuation.  We  find'  that  both  are  functions  of  tempera¬ 
ture  in  Figs.  9  through  2  2,  A  drastic  change  of  curvature  is  seen 
around  600°C,  which  is  half  the  melting  point  (60Q-650°C)  for  Remendur 
wire.  This  agrees  with  similar  data  for  other  alloys.  9 

* 

2.  Run  3,  Temperature  Increasing,  Modulus  Data 

For  Run  3  we  heated  the  samples  from  room  temperature  to 
920  C  and  drove  the  two  samples  with  magnetostrictive  transducers  of 
the  same  frequency  (140  kHz).  The  modulus  data  are  shown  in  Figs.  9 
and  10  with  square  symbols  for  the  data  points.  For  longitudinal  modulus , 
we  find  the  two  curves  are  of  the  same  shape  and  the  difference  in  value 
is  within  our  experimental  errors. 

3.  Run  4,  Temperature  Increasing,  Modulus  Data 

In  this  experiment,  we  cut  the  length  of  the  transducer  stub  for 
Sample  3  to  less  than  half  of  its  original  length  to  make  the  frequency 
of  the  driven  pulse  200  kHz,  but  for  Sample  4  we  keep  the  same  fre- 
r  '.y  (140  kHz).  Temperatures  are  from  room  temperature  to  927°C 
data  points  are  circles.  The  two  cuvvec  for  modulus  are  of  the 
same  shape  and  the  difference  is  within  experimental  errors. 

4.  Run  4,  Temperature  Decreasing,  hco^-o'us  Data 

The  measurements  are  taken  as  the  temperature  in  Run  4  returned 
to  room  temperature.  In  comparing  the  two  curves  for  Sample  4  and 
Sample  3,  we  find  3  percent  difference  between  the  two.  In  comparing 
Run  4  down  with  Run  4  up  for  each  sample,  we  find  that  for  "■'mple  4  they 
go  almost  the  same  curve  except  the  hysteresis,  in  the  range  . ..  900°C. 

However,  for  Sample  3  the  difference  between  Run  4  up  and  Run  4  down 
?.s  evident  and  difference  is  up  to  4  percent  at  temperatures  below  500°C. 
We  suspect  that  this  difference  is  caused  by  the  changing  of  the  length  of 
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the  transducer  stub  of  the  Remendur  for  Sample  3.  Further  experiments 
are  planned  to  study  this.  If  this  is  true,  the  longitudinal  modulus  is  a 
function  of  frequency  as  well  as  temperature. 

5.  Reflection  Coefficient 


For  reflection  coefficient,  we  found  it  to  be  around  60%  through 
all  the  experiments. 

6.  Attenuation 


In  general,  attenuation  began  to  increase  drastically  above  half 
the  melting  point  absolute.  See  Figs.  11  and  12.  At  lower  temperatures, 
we  find  many  negative  points  as  shown  in  Fig.  11,  but  when  we  calculate 
the  maximum  errors,  we  find  that  our  maximum  error  can  be  very  big 
with  respect  to  small  values  of  a .  (The  calculation  is  in  IIXC.  )  Reason¬ 
able  curves  can  be  drawn  if  the  maximum  errors  are  considered.  We 
found  that  a  is  the  same  (within  experimental  errors)  below  500°C  for 
all  the  runs  with  both  samples  except  for  Run  4  down  with  Sample  4.  For 
Sample  4,  the  curve  we  get  for  Run  4  down  is  strange;  we  do  not  know 
the  reason  yet. 

D.  Statistical  Errors  in  Attenuation  and  Reflection  Coefficient 

In  the  pulse-echo  technique  with  a  lead-in  line,  to  determine  reflec¬ 
tion  coefficient  and  attenuation  coefficient  we  use  the  amplitudes  of  echoes, 
the  interface  echo  A  and  two  free  back  echoes,  B  and  C.  Appropriate 
formula  has  been  derived  by  E.  P.  Papadakis*®’  ** 

R  =  [  AC/ (AC  -  1)1  1/2 

«  =  [inf  R  /C  )]  /2L 

where  A  =  —  B  =  ~ 

In  measuring  the  A,  B,  and  C  echoes,  our  errors  can  be  +  1%,  so,  when 
error  is  taken  into  consideration 

A  — >  A  (1  +  6  )  B  — >  B  (1  +  6  )  C - ->  (1  +  6  ) 

where  6  =  C.  01 
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I 


*  =FTT tty  =b.<1±26:>  =F<1ie>-  ■ 


° •§»**'> -io*'*"- 


where  e  =  0.  02 


r  .=4£- 

A  C  -  1 


=>  /  ^  (1  jgj>  C(l  +,e  )  1 

}a(1  +€)Ct(l  +e) -1  | 

(.  --  U/2 

=  L4giL±LiL . \ 

V-  C  (1  +  2  c  )  -  1, 


( _  )l/2 

J  A  C  (1  ±  2  6  )  V  , 

I  (A  C  -  1 )  <  1  +  44A-1 
l  [  -  AC  -  1  ) 


\1±Zt-u3^ 

C-l  L  +  A  C-l 


=  R  {  1  +  2  e  (l  -  R  )' 

2  { 

*  R  (  1  +  C  (1  -  R  }  1 

l  "  '! 

so  the  maximum  error  in  R  is  J 


max  f  _  .  2  v  -  _  _ _ 

5 -  =  e  (1  -  R  )  <  0.  02  -• 


The  statistical  deviation  is  calculated  for  a  as  follows: 


i 


•  ln  <-§■> 

a  =  -  the  R  and  C  are  different  in  sign 

2L 


2  sL  =  ln  (-  =) 

c 


2  ttL 


ln(e2“L)  =  ln(-~)  =  In  R  -  ln  C 

C 


2  a  L) 


dt  u  J 

_ (e_ _ 

2  aL 
e 


From  (18)  and  (20), 


2  «Lj 

2  aL 


%“l  +  l^pl  <  e  +  e  =0.04 


2  ^ 

A  (e  )  can  be  expressed  by  use  of  the  differentiation  formula 

d  ,  u.  u 
— —  (e  )  =  e 
du 

d  (eU)  =  eU  du 


A  (e U)  =  eU  Au 


A  ,  2  aL.  2  aL  . 

A  (e  )  =  e  A  (2  aL) 


so  Eq.  (21)  reduces  to 


A  (2  aL)  =  +  !~^|<2e  =  0.04 

R  C 


Since  L  is  held  constant, 


A  a  ** 


0.  02 


(23) 


Using  the  dimensions  in  Table  IV,  in  Sample  3  A  a  *  0.  003,  and  a  *  10 
at  temperature  below  500°C,  so  percentage  error  can  be  3000%.  That  is 
the  reason  that  a  comes  out  to  be  negative,  but  when  the  deviation  is  taken 
into  consideration  we  can  draw  a  reasonable  curve  with  a  *  0.  (There  is 
also  a  curve  of  percentage  error!*  in  R  and  2oL  with  respect  to  R  in 
Fig.  13.) 


IV.  PUBLICATIONS,  TALKS,  COURSES,  COMMITTEES,  STUDENTS 
A.  Papers  Published 


The  following  papers  on  work  supported  wholly  or  in  part  by  ONR 
have  been  published: 

E.  P.  Papadakis,  "On  the  Absolute  Accuracy  of  the  Pulse-Echo-Gverlap 
Method  and  the  Pulse-Superposition  Method  for  Ultrasonic  Velocity," 

J.  Acoust.  Soc.  Amer.  52,  843-846  (1972). 

E.  P.  Papadakis,  "Ultrasonic  Diffraction  Loss  and  Phase  Change  for 
Broadband  Pulses,"  J.  Acoust.  Soc.  Amer.  52,  847-S49  (1972). 

E.  P,  Papadakis,  L.  C.  Lynnworth,  K.  A.  Fowler,  and  E.  H.  Carnevale, 
"Ultrasonic  Attenuation  and  Velocity  in  Hot  Specimens  by  the  Momentary 
Contact  Method  With  Pressure  Coupling,  and  Some  Results  on  Steel  to 
1200C,"  J.  Acoust.  Soc.  Amer.  52,  850-857  (1972). 

E.  P.  Papadakis,  "Tabulation  of  the  Coefficients  of  a  Quadratic  Function 
for  the  Thermal  Expansion  of  Various  Alloys  and  Other  Engineering 
Materials,"  Materials  Science  and  Engineering  U),  195-203  (1972). 

E.  P.  Papadakis,  K.  A.  Fowler,  and  L.  C.  Lynnworth,  "New  Uses  for 
Ultrasonic  Spectrum  Analysis,  "  Proc.  1972  IEEE  Ultrasonics  Symposium, 
IEEE  #72  CHO  708-8  SU,  IEEE,  345  E  47  St. ,  N.  Y. ,  N.  Y.  10017,  1972, 

pp.  81-86. 

L.  C.  Lynnworth,  "Torsional  Wave  Transduction  in  Magentostrictive  Wire, 
Ultrasonics  10,  (5),  195-197  (Sept.  1972). 

L.  C.  Lynnworth,  "High  Temperature  Measurements,  "  Research/Develop- 
nnent  23[,  (10),  8  (Oct.  1972).  Note:  Over  100  responses  were  received 
concerning  this  Letter  to  the  Editor. 

B.  Papers  Scheduled 

The  following  papers  on  work  supported  wholly  or  in  part  by  ONR 
«.re  scheduled  for  publication  in  the  near  future: 

L.  C.  Lynnworth,  E.  P.  Papadakis,  and  K.  A.  Fowler,  "Ultrasound  Propa 
gation  Measurements  and  Applications,"  Proc.  ARPA-AMMRC  Symposium 
on  Advanced  NDT  Techniques  (to  be  published). 
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K.  A.  Fowler,  "Ultrasonic  Transducer  Characterization,'1  Proc,  ARPA- 
AMMRC  Symposium  on  Advanced  NDT  Techniques  (to  be  published). 

L.  C.  Lynnworth  and  E.  H.  Carnevale,  "Ultrasonic  Thermometry  Using 
Pulse-Echo  Techniques,"  Proc,  5th  Symposium  on  Temperature-Its 
Measurement  and  Control  in  Science  and  Industry  (to  be  published). 

L.  C.  Lynnworth,  "Attenuation  Measurements  Using  the  Pulse-Echo  AB 
Method,"  Materials  Evaluation  (to  be  published). 

E.  P.  Papadakis,  K.  A.  Fowler,  and  L.  C.  Lynnworth,  "Ultrasonic 
Attenuation  by  Spectrum  Analysis  of  Pulses  in  Buffer  Rods:  Method  and 
Diffraction  Corrections,"  J.  Acoust.  Soc.  Amer.  (to  be  published). 

E.  P.  Papadakis,  "Lens  Equation  for  Focused  Transducers, "  International 
Journal  of  Nondestructive  Testing  (to  be  published). 

E.  P.  Papadakis,  "Balanced  Resonator  for  Infrasonic  Measurements  of 
Young's  Modulus  and  Damping  in  Flexure,"  ASTM  J.  Testing  and  Evaluation 
(to  be  published). 

L.  C.  Lynnworth,  "Ultrasonic  Measurement  of  Elastic  Moduli  in  Slender 
Specimens  Using  Extensional  and  Torsional  Wave  Pulses,"  ASTM  J.  Testing 
and  Evaluation  (to  be  published). 

E.  P.  Papadakis,  "Measurement  of  Small  Changes  in  Ultrasonic  Velocity 
and  Attenuation,"  (Invited)  Crit.  Rev.  Solid  State  Sciences  (to  be  published). 

E.  P.  Papadakis,  "Ultrasonic  Diffraction  From  Single  Apertures,"  Physical 
Acoustics,  Principles  and  Methods,  W.  P.  Mason  and  R.  N.  Thurston,  Eds. 
Academic  Press,  N.  Y.  (future  volume  to  be  published). 

C.  Talks  Given 


The  following  talks  on  work  supported  wholly  or  in  part  by  ONR 
were  given  recently: 

E.  P.  Papadakis,  "Ultrasonic  Methods  for  Modulus  Measurement  in  Paper," 
TAPPI  Testing  Conference,  Philadelphia,  Pa.  ,  Sept.  26-29,  1972. 

E.  P.  Papadakis,  K.  A.  Fowler,  and  L.  C.  Lynnworth,  "New  Uses  for 
Ultrasonic  Spectrum  Analysis,"  IEEE  Ultrasonics  Symposium,  Boston, 

Mass. ,  Oct.  3-7,  1972. 
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C.  A.  Carey,  "Sound  Propagation  in  High-Temperature  Plasmas,"  IEEE 
Ultrasonics  Symposium,  Boston,  Mass. ,  Oct.  3-7,  1972. 


L.  C.  Lynnworth  and  N.  E.  Pedersen,  "Ultrasonic  Mass  Flowmeter, " 
IEEE  Ultrasonics  Symposium,  Boston,  Mass. ,  Oct.  3-7,  1972. 

D.  Talks  Scheduled 


The  following  talks  on  work  supported  wholly  or  in  part  by  ONR  are 
scheduled  for  the  near  future: 

L.  C.  Lynnworth,  "Ultrasonic  Temp.  Profiling  and  Attenuation, "  ASTM- 
E20. 06,  Washington,  D.  C.  (Nov.  17,  1972). 

L.  C.  Lynnworth,  "Nuclear  Reactor  Applications  of  Ultrasound, "  IEEE 
Boston  Section,  Joint  Meeting  of  Groups  on  Nuclear  Science  and  on  Sonics/ 
Ultrasonics  (Feb.  5,  1973). 

E.  Professional  Societies 


L.  C.  Lynnworth  is  secretary  of  ASTM  Subcommittee  E -20.  06, 
Acoustical  Thermometry.  He  i6  also  chairman  of  the  Boston  Section  of 
the  IEEE  Group  on  Sonics  and  Ultrasonics.  E.  P.  Papadakis  is  Editor  of 
the  Newsletter  and  Associate  Editor  for  industrial  ultrasonics  of  the  Trans¬ 
actions  on  Sonics  and  Ultrasonics  of  the  IEEE  Group  on  Sonics  and  Ultra¬ 
sonics  with  about  1200  members.  He  is  also  serving  on  the  Physical  Acoustics 
committee  of  the  Acoustical  Society  of  America,  on  the  Working  Group  on 
Ultrasonic  Attenuation  and  Velocity  (E-07.  06.  13)  of  ASTM  Subcommittee 
E-07.  06,  on  ASTM  Subcommittee  D-20.  15  concerned  with  Dynamic  Properties 
of  Plastics,  on  an  ad-hoc  committee  for  improved  nondestructive  testing 
methods  under  the  auspices  of  the  National  Science  Foundation  and  the  Air 
Force  Materials  Laboratory,  and  on  the  independent  Acoustic  Emission 
Working  Group.  K.  A.  Fowler  is  also  serving  on  the  latter  group. 

F.  Students 

James  Chen,  a  graduate  student  in  the  Physics  Dept,  at  Boston 
College,  is  working  at  Panametrics  on  this  research.  It  is  anticipated  that 
he  will  find  a  doctoral  thesis  topic  within  the  scope  of  the  work  under  this 
contract. 

G.  Course 

Our  course  #2911  Ultrasonics,  will  be  offered  again  at  the  North¬ 
eastern  University  Extension  School  when  sufficient  interest  has  been  shown 
to  warrant  it  in  the  judgment  of  the  University  administration. 
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Coil  Multiple  Segments 


Fig.  3.  I'ine  containing  four  sensing  segments:  Lj  =  L3  =  60  cm,  L2  =  L4  =  40  cm, 
Lead-in  and  segments  are  all  of  one  material.  Diameter  steps  are  designed  to 
provide  echoes  large  enough  at  room  temperature  to  compensate  for  attenuation 
at  high  temperature.  Lengths  L.  are  staggered  to  avoid  echo  interferences  at 
room  temperature. 
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Fig.  7.  1000  C  furnace  configuration. 
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Fig.  8.  Rack  configuration  for  holding  lead-in  wires  and 
thermocouple. 
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MONITORING  THE  MODULI  OF  POLYMERS 
WITH  ULTRASOUND 


by- 

Emmanuel  P,  Papadakis 


ABSTRACT 


Two  ultrasonic  methods  are  presented  for  the  measurement 
and  monitoring  of  the  elastic  moduli  of  polymers  and  other  liquids 
undergoing  solidification.  The  methods  are  thought  to  be  most  sensitive 
to  the  build-up  of  the  real  component  of  the  modulus  of  rigidity  of  the 
material.  Ultrasonic  traveling  waves  are  used  to  interrogate  a  boundary 
between  the  solid  body  transmitting  the  wave  and  the  sample  under  test. 
Data  on  epoxy  resins  are  given. 
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I.  INTRODUCTION 


The  behaviour  of  epoxy  resins  during  hardening  was  studied  as  a 
function  of  time  by  means  of  ultrasonics.  Both  bulk  wave  and  wire  wave 
methods  were  employed,  and  were  found  to  complement  each  other.  The 
propagation  media  (glass  blocks  and  tungsten  wires)  were  partially  immer¬ 
sed  in  beakers  of  epoxy,  and  the  increasing  energy  transfer  into  the  epoxy 
noted  as  it  hardened.  Temperature  was  monitored  simultaneously.  Quanti 
tative  results  were  obtained,  but  their  interpretation  is  complicated  by  two 
factors: 

1)  The  exothermic  nature  of  the  polymerization  and  cross-linking  reactions 
which  raised  the  temperature  in  the  epoxy,  and 

2)  The  fact  that  the  glass  transition  temperature  rises  during  polymeri¬ 
zation,  being  below  the  epoxy  temperature  at  the  start  of  polymerization 

(1.  2) 

and  above  the  epoxy  temperature  after  curing. 

« 

The  methods,  however,  hold  promise  for  further  study  of  polymers 
and  plastics. 

II.  EXPERIMENTS  WITH  BULK  WAVES 

A  glass  buffer  rod  with  plane  parallel  ends  having  an  ultrasonic 
transducer  bonded  to  one  end  was  inserted  into  a  pot  of  epoxy  mixed  with 
hardener  as  shown  in  Fig.  1.  Echoes  within  the  rod  were  observed  on  an 
oscilloscope  and  recorded  photographically  during  the  curing  process. 

The  temperature  of  the  epoxy  was  monitored  by  a  thermocouple,  and  is 
plotted  as  a  function  of  log  (time)  in  Fig.  2,  along  with  the  temperature 
of  the  curing  oven.  The  portion  of  the  wave  entering  the  epoxy  did  not 
return  to  the  sample  because  of  attenuation  in  the  epoxy  and  diffusion  at 


the  uneven  surface  of  the  pot.  Thus  the  amplitudes  of  the  first  and  second 
echoes  in  the  glass  sample  give  the  reflection  coefficient  between  the  glass 
and  the  epoxy  as 


R 


VA! 

<VAl'o 


(1) 


where  the  denominator  with  the  subscript  £  is  the  amplitude  ratio  before 
the  sample  is  dipped  into  the  pot  of  epoxy.  This  normalization  takes  care 
of  attenuation  and  beam  spreading  in  the  buffer  rod,,  R  is  related  to  the 

( 3) 

mechanical  impedances  of  the  glass  Z  and  epoxy  Z  by  the  formula 

G  & 
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Vz 

vz 


E 
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which  can  be  solved  for  the  unknown  Z_  as  follows: 

E 


(2) 


Z 


1-R 
G  1+R 


(3)* 


Z  is  found  as  a  function  of  time  by  recording  A  and  A  as  functions  of 
E  i  c 

time.  Since  density  p  remains  almost  constant  during  the  reaction,  and 

since^  Z  =  pv  where  v  is  the  ultrasonic  velocity,  Z  is  essentially  a 

measure  of  v.  Since^  v  =  (M/p)^2  where  M  is  the  appropriate  modulus 

1/2 

of  elasticity  for  the  wave  motion,  Z E  is  really  a  measure  of  (M)  as  it 
exists  in  the  epoxy  at  any  instant. 

DER  332  epoxy  resin'  and  MPDA  hardener  mixed  6.  5:1  by  weight 
were  used  in  this  experiment  with  5  MHz  ultrasonic  waves,  both  shear 
and  longitudinal.  The  cure  was  effected  in  an  oven  held  at  6l  +  1  C.  The 
data  on  the  reflection  coefficients  versus  time  are  plotted  in  Fig.  3.  In 
the  longitudinal  case,  the  reflection  coefficient  dropped  from  an  initial 
value  of  0.  77  at  the  glass-liquid  interface  to  a  final  value  of  0.  620  at  the 
glass-cured  solid  interface.  The  time  for  half  the  change  to  occur  was 
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200  min.  In  the  shear  wave  case,  R  dropped  from  an  initial  value  of 
0.  990  to  a  final  value  of  0.  715  with  a  half-life  of  170  min.  The  difference 
in  curing  times  is  due  to  differences  in  time  and  temperature  during  mix¬ 
ing  and  outgassing  the  mixtures  of  res:n  and  hardener  before  applying  them 
to  the  ultrasonic  samples. 

One  notices  that  the  change  in  reflectivity  occurs  in  a  precipitous 

and  symmetrical  manner  about  the  time  of  half-hardness  when  plotted  on 

semilogarithmic  paper  in  Fig.  3.  These  data  were  used  to  compute  ZE 

from  Eq.  (3).  The  relative  change  of  Z  with  time  is  shown  in  Fig.  4. 

E 

The  growth  is  very  rapid  about  t^^.  limiting  values  of  Z^,  are  given 

in  Table  I  in  terms  of  the  impedance  of  the  glass  found  from  measured 

3  5  5 

values  of  p  =  2.  21  g/cm  ,  v  =  5.  6l  x  10  cm/ sec,  and  v  =  3.  44  x  10 
'  I  s 

cm/sec. 


Table  I 

o  XX  2 

Impedances  of  DER  332  at  6l  C  in  10  gm  cm  sec 


Wave 

Mixed 

Cured 

Z_  /  Z 

E  G 

ZE 

VZG 

ZE 

Longitudinal 

0.  123 

1.  52 

0.  237 

2.  94 

Shear 

0.  005 

0.  038 

0.  166 

1.  26 

III.  EXPERIMENTS  WITH  WIRE  WAVES 


Three  experiments  were  performed  with  ultrasonic  waves  in  wires 
to  test  single-ended  wire  delay  lines  as  monitors  for  polymerization.  In 
these  cases,  the  wavelength  is  much  larger  than  the  wire  diameter,  so 
longitudinal  waves  travel  in  the  long  wavelength  extensional  wave  limit,  and 
torsional  waves  travel  in  a  single  mode.  All  three  experiments  utilized 
the  same  electronic  system,  a  Panatherm  50 10®  ultrasonic  thermometer 


instrument  with  a  magentostrictive  transducer,  a  wire  delay  line,  and  a 
monitoring  oscilloscope.  See  the  block  diagram  in  Fig.  5.  The  Panatherm 
pulses  the  transducer  coil  repetitively  with  a  current  spike,  generating 
stress  pulses  in  the  lead-in  wire.  These  pulses  echo  from  both  ends  of  a 
sensor  welded  to  the  end  of  the  lead-in  wire,  and  yield  echoes  seen  in 
Fig.  6.  The  Panatherm  measures  the  time  between  selected  pairs  of  echoes, 
and  the  amplitude  of  these  echoes.  The  amplitudes  may  also  be  photographed 
on  the  monitoring  oscilloscope,  as  done  in  these  experiments.  From  the 
travel  time  between  echoes,  and  the  length  of  the  sensor,  one  can  compute 
the  velocity  in  the  sensor  material.  With  a  calibrated  sensor  (travel  time 
versus  temperature)  the  instrument  is  an  ultrasonic  thermometer.  From  the 
amplitudes  of  the  first  three  echoes,  one  may  compute^  the  attenuation  in 
the  sensor  and  the  reflection  coefficient  at  the  lead-in/sensor  weld. 

In  one  experiment,  the  sensor  was  immersed  vertically  in  a  vial  of 
epoxy  such  that  the  meniscus  coincided  with  the  lead-in/sensor  weld.  The 
purpose  was  to  measure  the  attenuation  in  the  wire  as  the  epoxy  became 

mere  viscous  and  then  hardened,  and  so  to  monitor  changes  in  the  epoxy. 

* 

The  echoes  were  photographed  from  the  oscilloscope  and  the  temperature 

of  the  epoxy  was  monitored  with  a  thermocouple.  A  Joule-Wiedemann 
(7) 

transducer  was  employed,  giving  both  extensional  and  torsional  waves. 
Reflection  coefficients  and  attenuation  were  computed,  and  are  shown  in 
Fig.  7.  It  can  be  seen  that  the  rapid  rise  in  attenuation  coincides  with  the 
rapid  rise  in  temperature  due  to  the  exothermic  reaction,  and  that  in  a 
matter  ot  minutes  the  attenuation  becomes  too  large  to  measure.  The 
attenuation  can  arise  from  viscous  damping  as  the  viscosity  increases,  and 
from  radiation  as  the  real  part  of  the  modulus  in  the  epoxy  grows. 

In  a  second  experiment,  a  lead-in  wire  with  no  sensor  attached  was 
immersed  to  a  length  of  about  10  cm  in  a  vial  of  epoxy  which  was  cooled  with 
tap  water  (8°C)  in  a  surrounding  coil.  The  purpose  was  to  observe  the  echo 
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which  was  found  to  be  generated  at  the  meniscus  when  the  epoxy  begins 

to  solidify.  On  solidification,  the  wire  within  the  epoxy  increases  in 

(4) 

effective  acoustic  impedance  because  of  the  adhering  epoxy,  and  a 
definite  echo  is  created  at  the  meniscus.  The  amplitude  of  this  echo 
was  followed  with  oscilloscope  photographs  as  the  epoxy  hardened.  A 
plot  is  shown  in  Fig.  8,  The  echo  becomes  visible  in  less  than  45  minutes, 
and  grows  slowly  at  this  temperature  until  a  rapid  increase  begins  at 
about  200  minutes.  The  climb  is  interrupted  at  point  A  on  the  graph  by 
an  apparent  partial  disbonding  between  the  wire  and  the  epoxy.  (Cracks 
were  visible  upon  inspection  after  the  experiment.  )  However,  the  echo 
amplitude  continues  to  climb  thereafter  until  another  discontinuity  appears 
at  point  B  on  the  graph.  It  is  thought  that  this,  too,  arose  from  partial 
disbonding.  The  echo  at  the  meniscus  is  thought  to  be  generated  by  the 
growth  of  the  shear  modulus  (modulus  of  rigidity)  of  the  epoxy. 

A  third  experiment  was  performed  to  test  the  automatic  capabilities 
of  the  Panatherm  5010  instrument  to  detect  and  record  the  initiation  and 
growth  of  the  echo  at  the  meniscus  of  the  epoxy.  The  system  was  set  up 
as  in  Fig.  9  with  the  wire  immersed  to  a  depth  2L  into  a  vial  of  epoxy. 

The  first  segment  of  length  L  was  a  sensor,  while  the  second  segment 
was  part  of  the  lead-in  wire.  The  "Selector  Monitor"  of  the  Panatherm 
was  set  to  observe  the  echoes  from  the  weld  (point  II)  and  the  end  of  the 
sensor  (point  III)  while  the  echo  from  the  meniscus  was  absent  or  small. 

The  amplitudes  A  and  B  as  measured  by  the  Panatherm  were  sent  to  a 
multi-channel  recorder.  The  trigger  level  in  the  Panatherm  was  set  to 
cause  the  instrument  to  begin  monitoring  the  echo  at  the  meniscus  (point  I 
of  Fig.  9)  when  it  had  grown  to  a  small  but  finite  value.  The  experiment 
when  run  resulted  in  data  sketched  schematically  in  Fig.  10  which  represents 
the  record  on  a  strip  chart.  Initially  the  echoes  from  points  II  and  III 
(charts  A  and  B)  fall  as  attenuation  increases.  When  the  echo  from  point  I, 


the  meniscus,  reaches  the  trigger  level,  the  instrument  jumps  to  monitor 
points  I  and  II  on  charts  A  and  B  instead.  The  meniscus  echo  grows  while 
the  echo  at  point  II  continues  to  decrease  in  amplitude. 

IV.  DISCUSSION 

Two  ultrasonic  methods  have  been  presented  for  studying  the  growth 
of  the  elastic  moduli  of  epoxy  resins  during  polymerization.  By  extension, 
the  methods  could  be  used  for  other  polymers  and  for  liquids  which  change 
to  a  solid  phase  at  a  melting  point  or  go  through  a  glass  transition.  The 
two  methods,  using  wire  waves  and  bulk  waves  respectively,  can  operate 
from  a  few  kHz  to  several  GHz.  The  dividing  line  between  them  is  in  the 
range  500  kHz  to  2  MHz,  since  bulk  wave  transducers  for  plane  waves 
become  excessively  large  below  500kHz  and  since  wire-sonic  transmission 
lines  become  difficult  to  build  above  2  MHz.  The  methods  are  applicable 
to  the  study  of  bulk  quantities  of  liquids  undergoing  solidification.  For 
thin  layers,  surface  waves,  Lamb  waves  and  interface  waves  might  be* 
preferable,  although  the  wire  waves  should  be  as  sensitive.  The  methods 
react  to  relatively  gross  changes  in  system  properties.  For  measuring 
smaller  changes,  see  Ref.  8. 
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CAPTIONS 


Fig.  1.  Block  diagram  of  the  experiment  using  ultrasonic  bulk  waves 
impinging  upon  epoxy  resin  at  normal  incidence. 

Fig.  2.  Temperature  in  tae  oven  and  in  the  epoxy  in  experiment  performed 
as  in  Fig.  1. 

Fig.  3.  Reflection  coefficients  at  the  buffer/epoxy  interface  in  two  experi¬ 
ments  performed  as  in  Fig.  1  on  curing  epoxy. 

Fig.  4.  Inpedance  ratios  computed  from  the  data  of  Fig.  3.  See  Table  I 
for  absolute  values. 

Fig.  5.  Block  diagram  of  the  instrumentation  used  for  waves  in  wires. 

Specimens  and  thermocouples  were  dipped  vertically  into  vials  of 
epoxy. 

Fig.  6.  Echoes  (upper)  and  Selector  Monitor  pips  (lower)  photographed  on 

the  CRO,  and  corresponding  specimen  geometry,  in  one  experiment 
performed  with  equipment  in  Fig.  5  with  specimen  immersed. 

Fig.  7.  Data  on  a  wire  specimen  in  epox,  during  curing. 

Fig.  8.  Amplitude  of  the  meniscus  echo  from  a  lead-in  wire  (without 
specimen)  held  in  epoxy  during  curing. 

Fig.  9.  Test  set  to  monitor  polymerization  automatically.  The  amplitude 
outputs  A  and  B  were  from  the  weld  (II)  and  the  specimen  end  (III) 
respectively  at  first.  When  the  men;scus  echo  (I)  grew  to  a  preset 
trigger  level  channel  A  began  monitoring  the  meniscus  (I)  and 
channel  B  switched  to  the  weld  (II). 

Fig.  10.  Data  from  the  test,  set  in  Fig.  9.  When  the  meniscus  echo  reached 
the  trigger  level,  channels  A  and  B  switched  from  positions  II  and 
III  to  positions  I  and  II  of  Fig.  9. 
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APPENDIX  B 

Attenuation  and  Reflection  Coefficient  Nomogram 

by 

Lawrence  C.  Lynnworth 


In  determining  the  attenuation  coefficient  a  (Np/unit  length)  using 
the  pulse-echo  buffer  rod  technique,  calculations  usually  involve  three 
echoes,  such  as  the  interface  echo  A  and  the  free  back  surface  echoes  B 
and  C.  Appropriate  formulas  have  been  derived  by  Papadakis/1^  taking 
into  account  the  sound  pressure  reflection  coefficient  R  at  the  buffer/ 
specimen  interface.  These  formulas  usually  are  written  in  terms  of 
normalized  echo  amplitude  ratios  such  as  A  =  A/B  and  C  =  C/B  with 
signs  preserved,  and  the  specimen  length  L: 

|  R|  =  [  A  C/(A  C  -  I)]  1/2  (1) 

and  _ 

a  =  [In  (R/C)3/2L  .  (2) 

The  sign  of  R  is  the  same  as  that  of  C,  in  the  present  notation.  As  a 
cross  check,  the  sign  of  R  can  be  found  from  the  impedances  Z^  and  Z^ 
of  the  buffer  and  specimen  as 

R  positive  when  Z_  >  Z,  ; 

2  1  (3) 

R  negative  when  Z>2<  ^1  * 

For  highest  accuracy,  A,  B  and  C  may  be  corrected  for  diffraction 

(2) 

prior  to  evaluating  a  and  R  by  these  formulas.  They  should  also  be 
corrected  for  nonlinearities  in  the  transducer  and  electronic  amplifiers. 

In  some  cases  it  is  convenient,  and  sufficiently  accurate,  to  de¬ 
termine  2a  L  and  R  graphically  from  a  log-log  nomogram.  One  such 
nomogram  is  constructed  by  plotting  iso-ar  curves  and  iso-R  lines  as 
functions  of  A  and  C.  As  a  numerical  example,  if  A  =  10  and  C  =  -0. 1, 
we  estimate  2aL  *  2  Np  and  |  R|  ~  0.  7.  (More  exact  values,  computed 
using  (1)  and  (2),  are:  |  Rj  =  0.  7071,  2a L  =  1.  956  Np  3  17  dB.  ) 
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Fig.  1.  Nomogram  for  finding  specimen  loss  and  reflection  coefficient 
at  the  end  of  a  buffer  rod  from  the  amplitudes  of  the  first  three  echoes 
as  corrected  for  beam  spreading  when  present,  and  transducer  and 
electronic  nonlinearity  when  present. 
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I  INTRODUCTION 


My  subject  today  is  the  modulus  of  paper,  and  two  ultrasonic 
methods  to  measure  it  on-lin^. 

Modulus  is  related  to  strength  and  to  bending  stiffness.  (M. 
Jackson  and  G.  Gavelin,  Svensk;  Pappers  Tidnjng  70,  63  (1970).) 

The  latter  is  proportional  to  Young' s  modulus,  and  the  former  can 
be  correlated  with  it.  It  would  be  advantageous  to  be  able  to  measure, 
on-line,  a  quantity  which  is  a  function  of  Young' s  modulus. 

II.  PLATE  WAVE  THEORY 

One  such  quantity  is  the  velocity  Vp  of  the  first  longitudinal 
plate  wave  in  the  paper.  This  wave  velocity  occurs  when  the  paper 
is  much  thinner  than  one  ultrasonic  wavelength. 

The  modes  in  plates  are  illustrated  in  Slide  1.  (First  slide, 
please. )  The  wave  velocity  is  in  the  direction  of  the  vector  V,  while 
the  particle  motion  directions  are  shown  by  the  short  arrows.  The 
first  longitudinal  mode,  the  first  flexural  mode,  and  thp  zeroth  shear 
mode  are  shown. 

The  group  delay  (or  travel  time)  of  modes  in  plates  are 
illustrated  in  Slide  2.  (Next  slide,  please. )  The  group  delay  D 
of  wave  pulses  is  plotted  against  frequency  f.  Both  quantities  are 
normalized  with  respect  to  thickness  and  velocity  to  be  dimensionless. 
The  plate  velocity  Vp  is  indicated  at  the  lower  end  of  the  first  longi¬ 
tudinal  mode  curve.  This  i3  the  fastest  mode  in  the  plate,  showing  the 
shortest  delay  D.  The  velocity  is  given  by  the  equation 

Vp  =  [  Y/p  (1  -  a2)]  1/2  (1) 

where  p  is  density,  or  is  Poisson' s  ratio,  and  Y  is  Young' s  modulus. 
Other  relevant  equations  are  shown  in  Slide  3.  (Next  slide,  please. ) 

This  formula  for  Vp  is  not  ideal  ’/ec<.use  it  depends  on  Poisson' s 
ratio  as  well  as  on  Young' s  modulus.  However,  Poisson' s  ratio  may 
be  expressed  as 


0  =  Y/2p  -  l  (2) 

\ 

where  ji  is  the  modulus  of  rigidity.  The  rigidity  modulus  may  be  found 
by  propagating  a  shear  wave  in  the  plate.  The  velocity  Vg  of  the  zeroth 
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Slide  1.  Illustration  of  motion 
of  lowest  modes  in  plates. 


Slide  2.  Group  delay  of  some 
low  modes  in  isotropic  plates. 


Slide  3.  Slide  3  contained 
Equations  (1)  through  (5)  in] 
the  text. 


Slide  4.  Wheel  transducers  with 
piezoelectric  elements  for  gen¬ 
erating  longitudinal  and  shear 
wave  beams  in  plates. 


Slide  5.  Instrumentation  for  velocity 
measurements  in  paper  with  wheel- 
type  transducers. 
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Slide  6.  Velocity  in  paper  for 
first  longitudinal  and  zeroth 
shear  modes  by  wheel  method. 


Slide  7.  Velocity  in  card  stock 
for  first  longitudinal  and  zeroth 
shear  modes  by  wheel  method. 


Slide  8.  Synchronous  wheel 
transducer  configuration  for 
measuring  velocity  in  sheet 
material  on-line. 


Slide  9.  Noncontacting  experiment 
on  velocity  in  paper.  Transmission 
through  the  paper  occurs  when  plate 
waves  are  excited  in  it.  (Luukkala 
et  al.  ,  Ultrasonics  9.,  201-208  (1971).) 


Slide  10.  Wavelength  matching 
condition  for  transmission 
through  the  paper.  (Luukkala 
et  al.  ,  op.  cit.  ) 


Slide  12.  Amplitude  versus  angle 
for  several  frequencies  by  the  non¬ 
contacting  method  in  card  0.  85  mm 
thick.  Maximum  signal  occurs  at 
the  wavelength  matching  condition 
which  excites  a  plate  wave. 
(Luukkala  et  al.  ,  op.  cit. ) 


Slide  11.  Slide  11  contained 
Equation  (8)  in  the  text. 


IOO  200  300 

Frequency  f  [kHz) 


Slide  13.  Mode  curves  of  phase 
velocity  versus  frequency  com¬ 
puted  by  Equation  (6)  from  the 
data  of  Slide  12.  (Luukkala  et  al.  , 
op.  cit. ) 


Slide  14.  Extrapolation  of  data  to 
low-frequency  limit  to  get  plate 
velocity  Vp  of  Equation  (1). 
(Luukkala  et  al. ,  op.  cit. ) 


Table  1  Comparison  between  the  two  methods  (0  =  machine  direction,  90  =  transverse  direction,  x  =  measurement  could  not 
be  carried  out) 


Plate  wave  resonance  Dynamic  modulus  tester 


Sample  code 

llo  , 
[ms"1] 

^Lgo, 

(ms-1) 

Al 

^L0  , 

[ms"1  ] 

'■'L90, 
(ms-1 ) 

Al 

SI-1 

1890 

1530 

1,23 

1980 

1570 

SI-2 

1930 

2070 

1720 

1,22 

SI -3 

1930 

1,26 

1830 

S2-1 

2200 

1490 

1,48 

3050 

2300 

1,33 

S5-1 

4800 

X 

X 

3530 

2020 

1,75 

S5-2 

3320 

X 

X 

1,9b 

S5-3 

2790 

X 

X 

3530 

1950 

1,81 

S6-1 

1710 

1.24 

2120 

1580 

1,34 

S6-2 

1960 

1,53 

2660 

1650 

1,61 

S6-3 

1720 

8jB6s  |TSTi»BiB 

1,15 

1610 

1540 

C185 

2760 

1870 

1,48 

3260 

2220 

1,47 

C74 

3020 

2000 

1,5 

3140 

2600 

1,21 

C427 

2670 

X 

X 

2980 

2490 

B638 

2780 

2410 

1,15 

3290 

2840 

1.16 

B1572 

3940 

2600 

1,51 

X 

X 

X 

K 

2790 

1,44 

4450 

2440 

1,82 

B1 

X 

X 

1380 

2,35 

B2 

2410 

1,70 

4200 

1820 

2,31 

B3 

3830 

2020 

1,91 

3770 

1940 

1,94 

Slide  15.  Extrapolated  results  compared  with  data  from 
dynamic  modulus  tester.  (Luukkala  et  al.  ,  op.  cit.  ) 


Slide  16.  Conceptual  sketch  of 
on-line  modulus  tester  for  sheet 
material  employing  noncontacting 
method. 


TRANSMITTING 


Slide  17.  Noncontacting  test  set-up 
with  two  pairs  of  transducers  at  two 
angles  and  radiating  at  two  frequen¬ 
cies  for  a  dispersive  wave  mode. 


Slide  18.  Contact  wheel  transducer 
radiating  in  untouched  paper  to  a 
noncontacting  receiver. 
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be  done  directly  by  rotating  the  transducers,  or  it  could  be  done  in¬ 
directly  by  varying  the  frequency  fed  to  a  phased  array  transducer. 
Physical  angle  or  frequency  would  have  to  be  varied  with  a  feedback 
mechanism  to  hunt  for  maximum  transmission.  The  hunt  for  maxi¬ 
mum  transmission  could  be  complicated  by  various  ambient  factors 
such  as  turbulent  air  or  absorption  in  the  paper. 

Ordinary  plane  transducers  could  be  used  in  the  noncontacting 
case  if  one  chose  to  operate  in  a  dispersive  region  of  the  longitudinal 
mode.  For  each  direction  of  interest,  e.  g.  ,  the  machine  direction, 
two  pairs  of  transducers  would  be  needed.  See  Slide  17.  (Next  slide, 
please.  )  The  two  pairs  would  be  set  at  different  angles  of  incidence, 
and  the  input  frequency  to  each  pair  would  be  different  and  would  be 
swept  independently  to  find  maximum  transmission  for  each  pair.  The 
resultant  velocity  of  plate  waves  would  be  found  by  extrapolating  a  plot 
of  f  sin^g  versus  f  to  zero  frequency  as  show-  .arlier  in  Slide  14.  A 
single  system  to  monitor  frequency  at  a  single  angle  might  be  adequate 
for  practical  purposes. 

The  contact  method  seems  attractive  at  first  glance  if  contact 
were  permissible.  Turbulent  air  problems  are  absent,  and  absorption 
in  the  paper  would  not  be  a  problem  as  long  as  signal  strength  did  not 
fall  below  a  threshold.  However,  the  rolling  action  under  pressure 
might  modify  the  paper  being  measured.  This  problem  might  be  over¬ 
come  with  an  air-coupled  transducer  as  receiver  "upstream"  from  the 
wheel  transmitter.  See  Slide  18.  (Next  slide,  please. )  The  virtue  of 
this  configuration  would  be  that  paper  untouched  by  the  wheel  roller 
would  be  interrogated. 

In  summary,  on-line  modulus  monitoring  of  paper  is  now  within 
the  realm  of  possibility.  Work  is  still  necessary  to  define  the  best 
method  to  use,  and  the  best  configuration  to  exploit  the  method.  A 
close  interaction  between  producers  and  experimenters  will  be  neces¬ 
sary  to  create  a  viable  system. 
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